Neural network emulation of the human ventricular cardiomyocyte action potential for more efficient computations in pharmacological studies.

Computer models of the human ventricular cardiomyocyte action potential (AP) have reached a level of detail and maturity that has led to an increasing number of applications in the pharmaceutical sector. However, interfacing the models with experimental data can become a significant computational burden. To mitigate the computational burden, the present study introduces a neural network (NN) that emulates the AP for given maximum conductances of selected ion channels, pumps, and exchangers. Its applicability in pharmacological studies was tested on synthetic and experimental data. The NN emulator potentially enables massive speed-ups compared to regular simulations and the forward problem (find drugged AP for pharmacological parameters defined as scaling factors of control maximum conductances) on synthetic data could be solved with average root-mean-square errors (RMSE) of 0.47 mV in normal APs and of 14.5 mV in abnormal APs exhibiting early afterdepolarizations (72.5% of the emulated APs were alining with the abnormality, and the substantial majority of the remaining APs demonstrated pronounced proximity). This demonstrates not only very fast and mostly very accurate AP emulations but also the capability of accounting for discontinuities, a major advantage over existing emulation strategies. Furthermore, the inverse problem (find pharmacological parameters for control and drugged APs through optimization) on synthetic data could be solved with high accuracy shown by a maximum RMSE of 0.22 in the estimated pharmacological parameters. However, notable mismatches were observed between pharmacological parameters estimated from experimental data and distributions obtained from the Comprehensive in vitro Proarrhythmia Assay initiative. This reveals larger inaccuracies which can be attributed particularly to the fact that small tissue preparations were studied while the emulator was trained on single cardiomyocyte data. Overall, our study highlights the potential of NN emulators as powerful tool for an increased efficiency in future quantitative systems pharmacology studies.

Digital Twinning of Cardiac Electrophysiology Models From the Surface ECG: A Geodesic Backpropagation Approach.

The eikonal equation has become an indispensable tool for modeling cardiac electrical activation accurately and efficiently. In principle, by matching clinically recorded and eikonal-based electrocardiograms (ECGs), it is possible to build patient-specific models of cardiac electrophysiology in a purely non-invasive manner. Nonetheless, the fitting procedure remains a challenging task. The present study introduces a novel method, Geodesic-BP, to solve the inverse eikonal problem. Geodesic-BP is well-suited for GPU-accelerated machine learning frameworks, allowing us to optimize the parameters of the eikonal equation to reproduce a given ECG. We show that Geodesic-BP can reconstruct a simulated cardiac activation with high accuracy in a synthetic test case, even in the presence of modeling inaccuracies. Furthermore, we apply our algorithm to a publicly available dataset of a biventricular rabbit model, with promising results. Given the future shift towards personalized medicine, Geodesic-BP has the potential to help in future functionalizations of cardiac models meeting clinical time constraints while maintaining the physiological accuracy of state-of-the-art cardiac models.

Neural network emulation of the human ventricular cardiomyocyte action potential: a tool for more efficient computation in pharmacological studies.

Computer models of the human ventricular cardiomyocyte action potential (AP) have reached a level of detail and maturity that has led to an increasing number of applications in the pharmaceutical sector. However, interfacing the models with experimental data can become a significant computational burden. To mitigate the computational burden, the present study introduces a neural network (NN) that emulates the AP for given maximum conductances of selected ion channels, pumps, and exchangers. Its applicability in pharmacological studies was tested on synthetic and experimental data. The NN emulator potentially enables massive speed-ups compared to regular simulations and the forward problem (find drugged AP for pharmacological parameters defined as scaling factors of control maximum conductances) on synthetic data could be solved with average root-mean-square errors (RMSE) of 0.47mV in normal APs and of 14.5mV in abnormal APs exhibiting early afterdepolarizations (72.5% of the emulated APs were alining with the abnormality, and the substantial majority of the remaining APs demonstrated pronounced proximity). This demonstrates not only very fast and mostly very accurate AP emulations but also the capability of accounting for discontinuities, a major advantage over existing emulation strategies. Furthermore, the inverse problem (find pharmacological parameters for control and drugged APs through optimization) on synthetic data could be solved with high accuracy shown by a maximum RMSE of 0.21 in the estimated pharmacological parameters. However, notable mismatches were observed between pharmacological parameters estimated from experimental data and distributions obtained from the Comprehensive in vitro Proarrhythmia Assay initiative. This reveals larger inaccuracies which can be attributed particularly to the fact that small tissue preparations were studied while the emulator was trained on single cardiomyocyte data. Overall, our study highlights the potential of NN emulators as powerful tool for an increased efficiency in future quantitative systems pharmacology studies.

A personalized real-time virtual model of whole heart electrophysiology.

Computer models capable of representing the intrinsic personal electrophysiology (EP) of the heart in silico are termed virtual heart technologies. When anatomy and EP are tailored to individual patients within the model, such technologies are promising clinical and industrial tools. Regardless of their vast potential, few virtual technologies simulating the entire organ-scale EP of all four-chambers of the heart have been reported and widespread clinical use is limited due to high computational costs and difficulty in validation. We thus report on the development of a novel virtual technology representing the electrophysiology of all four-chambers of the heart aiming to overcome these limitations. In our previous work, a model of ventricular EP embedded in a torso was constructed from clinical magnetic resonance image (MRI) data and personalized according to the measured 12 lead electrocardiogram (ECG) of a single subject under normal sinus rhythm. This model is then expanded upon to include whole heart EP and a detailed representation of the His-Purkinje system (HPS). To test the capacities of the personalized virtual heart technology to replicate standard clinical morphological ECG features under such conditions, bundle branch blocks within both the right and the left ventricles under two different conduction velocity settings are modeled alongside sinus rhythm. To ensure clinical viability, model generation was completely automated and simulations were performed using an efficient real-time cardiac EP simulator. Close correspondence between the measured and simulated 12 lead ECG was observed under normal sinus conditions and all simulated bundle branch blocks manifested relevant clinical morphological features.

GEASI: Geodesic-based earliest activation sites identification in cardiac models.

The identification of the initial ventricular activation sequence is a critical step for the correct personalization of patient-specific cardiac models. In healthy conditions, the Purkinje network is the main source of the electrical activation, but under pathological conditions the so-called earliest activation sites (EASs) are possibly sparser and more localized. Yet, their number, location and timing may not be easily inferred from remote recordings, such as the epicardial activation or the 12-lead electrocardiogram (ECG), due to the underlying complexity of the model. In this work, we introduce GEASI (Geodesic-based Earliest Activation Sites Identification) as a novel approach to simultaneously identify all EASs. To this end, we start from the anisotropic eikonal equation modeling cardiac electrical activation and exploit its Hamilton-Jacobi formulation to minimize a given objective function, for example, the quadratic mismatch to given activation measurements. This versatile approach can be extended to estimate the number of activation sites by means of the topological gradient, or fitting a given ECG. We conducted various experiments in 2D and 3D for in-silico models and an in-vivo intracardiac recording collected from a patient undergoing cardiac resynchronization therapy. The results demonstrate the clinical applicability of GEASI for potential future personalized models and clinical intervention.

A Framework for the generation of digital twins of cardiac electrophysiology from clinical 12-leads ECGs.

Cardiac digital twins (Cardiac Digital Twin (CDT)s) of human electrophysiology (Electrophysiology (EP)) are digital replicas of patient hearts derived from clinical data that match like-for-like all available clinical observations. Due to their inherent predictive potential, CDTs show high promise as a complementary modality aiding in clinical decision making and also in the cost-effective, safe and ethical testing of novel EP device therapies. However, current workflows for both the anatomical and functional twinning phases within CDT generation, referring to the inference of model anatomy and parameters from clinical data, are not sufficiently efficient, robust and accurate for advanced clinical and industrial applications. Our study addresses three primary limitations impeding the routine generation of high-fidelity CDTs by introducing; a comprehensive parameter vector encapsulating all factors relating to the ventricular EP; an abstract reference frame within the model allowing the unattended manipulation of model parameter fields; a novel fast-forward electrocardiogram (Electrocardiogram (ECG)) model for efficient and bio-physically-detailed simulation required for parameter inference. A novel workflow for the generation of CDTs is then introduced as an initial proof of concept. Anatomical twinning was performed within a reasonable time compatible with clinical workflows (<4h) for 12 subjects from clinically-attained magnetic resonance images. After assessment of the underlying fast forward ECG model against a gold standard bidomain ECG model, functional twinning of optimal parameters according to a clinically-attained 12 lead ECG was then performed using a forward Saltelli sampling approach for a single subject. The achieved results in terms of efficiency and fidelity demonstrate that our workflow is well-suited and viable for generating biophysically-detailed CDTs at scale.

Automatic reconstruction of the left atrium activation from sparse intracardiac contact recordings by inverse estimate of fibre structure and anisotropic conduction in a patient-specific model.

AIMS: Electric conduction in the atria is direction-dependent, being faster in fibre direction, and possibly heterogeneous due to structural remodelling. Intracardiac recordings of atrial activation may convey such information, but only with high-quality data. The aim of this study was to apply a patient-specific approach to enable such assessment even when data are scarce, noisy, and incomplete. METHODS AND RESULTS: Contact intracardiac recordings in the left atrium from nine patients who underwent ablation therapy were collected before pulmonary veins isolation and retrospectively included in the study. The Personalized Inverse Eikonal Model from cardiac Electro-Anatomical Maps (PIEMAP), previously developed, has been used to reconstruct the conductivity tensor from sparse recordings of the activation. Regional fibre direction and conduction velocity were estimated from the fitted conductivity tensor and extensively cross-validated by clustered and sparse data removal. Electrical conductivity was successfully reconstructed in all patients. Cross-validation with respect to the measurements was excellent in seven patients (Pearson correlation r > 0.93) and modest in two patients (r = 0.62 and r = 0.74). Bland-Altman analysis showed a neglectable bias with respect to the measurements and the limit-of-agreement at -22.2 and 23.0 ms. Conduction velocity in the fibre direction was 82 ± 25 cm/s, whereas cross-fibre velocity was 46 ± 7 cm/s. Anisotropic ratio was 1.91±0.16. No significant inter-patient variability was observed. Personalized Inverse Eikonal model from cardiac Electro-Anatomical Maps correctly predicted activation times in late regions in all patients (r = 0.88) and was robust to a sparser dataset (r = 0.95). CONCLUSION: Personalized Inverse Eikonal model from cardiac Electro-Anatomical Maps offers a novel approach to extrapolate the activation in unmapped regions and to assess conduction properties of the atria. It could be seamlessly integrated into existing electro-anatomic mapping systems. Personalized Inverse Eikonal model from cardiac Electro-Anatomical Maps also enables personalization of cardiac electrophysiology models.

Learning atrial fiber orientations and conductivity tensors from intracardiac maps using physics-informed neural networks.

Electroanatomical maps are a key tool in the diagnosis and treatment of atrial fibrillation. Current approaches focus on the activation times recorded. However, more information can be extracted from the available data. The fibers in cardiac tissue conduct the electrical wave faster, and their direction could be inferred from activation times. In this work, we employ a recently developed approach, called physics informed neural networks, to learn the fiber orientations from electroanatomical maps, taking into account the physics of the electrical wave propagation. In particular, we train the neural network to weakly satisfy the anisotropic eikonal equation and to predict the measured activation times. We use a local basis for the anisotropic conductivity tensor, which encodes the fiber orientation. The methodology is tested both in a synthetic example and for patient data. Our approach shows good agreement in both cases and it outperforms a state of the art method in the patient data. The results show a first step towards learning the fiber orientations from electroanatomical maps with physics-informed neural networks.

An Inverse Eikonal Method for Identifying Ventricular Activation Sequences from Epicardial Activation Maps.

A key mechanism controlling cardiac function is the electrical activation sequence of the heart's main pumping chambers termed the ventricles. As such, personalization of the ventricular activation sequences is of pivotal importance for the clinical utility of computational models of cardiac electrophysiology. However, a direct observation of the activation sequence throughout the ventricular volume is virtually impossible. In this study, we report on a novel method for identification of activation sequences from activation maps measured at the outer surface of the heart termed the epicardium. Conceptually, the method attempts to identify the key factors governing the ventricular activation sequence - the timing of earliest activation sites (EAS) and the velocity tensor field within the ventricular walls - from sparse and noisy activation maps sampled from the epicardial surface and fits an Eikonal model to the observations. Regularization methods are first investigated to overcome the severe ill-posedness of the inverse problem in a simplified 2D example. These methods are then employed in an anatomically accurate biventricular model with two realistic activation models of varying complexity - a simplified trifascicular model (3F) and a topologically realistic model of the His-Purkinje system (HPS). Using epicardial activation maps at full resolution, we first demonstrate that reconstructing the volumetric activation sequence is, in principle, feasible under the assumption of known location of EAS and later evaluate robustness of the method against noise and reduced spatial resolution of observations. Our results suggest that the FIMIN algorithm is able to robustly recover the full 3D activation sequence using epicardial activation maps at a spatial resolution achievable with current mapping systems and in the presence of noise. Comparing the accuracy achieved in the reconstructed activation maps with clinical data uncertainties suggests that the FIMIN method may be suitable for the patient- specific parameterization of activation models.

Correction to: A systematic literature review of simulation models for non-technical skill training in healthcare logistics.

[This corrects the article DOI: 10.1186/s41077-018-0072-7.].

A systematic literature review of simulation models for non-technical skill training in healthcare logistics.

Background: Resource allocation in patient care relies heavily on individual judgements of healthcare professionals. Such professionals perform coordinating functions by managing the timing and execution of a multitude of care processes for multiple patients. Based on advances in simulation, new technologies that could be used for establishing realistic representations have been developed. These simulations can be used to facilitate understanding of various situations, coordination training and education in logistics, decision-making processes, and design aspects of the healthcare system. However, no study in the literature has synthesized the types of simulations models available for non-technical skills training and coordination of care. Methods: A systematic literature review, following the PRISMA guidelines, was performed to identify simulation models that could be used for training individuals in operative logistical coordination that occurs on a daily basis. This article reviewed papers of simulation in healthcare logistics presented in the Web of Science Core Collections, ACM digital library, and JSTOR databases. We conducted a screening process to gather relevant papers as the knowledge foundation of our literature study. The screening process involved a query-based identification of papers and an assessment of relevance and quality. Results: Two hundred ninety-four papers met the inclusion criteria. The review showed that different types of simulation models can be used for constructing scenarios for addressing different types of problems, primarily for training and education sessions. The papers identified were classified according to their utilized paradigm and focus areas. (1) Discrete-event simulation in single-category and single-unit scenarios formed the most dominant approach to developing healthcare simulations and dominated all other categories by a large margin. (2) As we approached a systems perspective (cross-departmental and cross-institutional), discrete-event simulation became less popular and is complemented by system dynamics or hybrid modeling. (3) Agent-based simulations and participatory simulations have increased in absolute terms, but the share of these modeling techniques among all simulations in this field remains low. Conclusions: An extensive study analyzing the literature on simulation in healthcare logistics indicates a growth in the number of examples demonstrating how simulation can be used in healthcare settings. Results show that the majority of studies create situations in which non-technical skills of managers, coordinators, and decision makers can be trained. However, more system-level and complex system-based approaches are limited and use methods other than discrete-event simulation.